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Abstract 

We introduce a combined macroscopic/microscopic transport approach em- 
ploying relativistic hydrodynamics for the early, dense, deconfined stage of the 
reaction and a microscopic non-equilibrium model for the later hadronic stage 
where the equilibrium assumptions are not valid anymore. Within this ap- 
proach we study the dynamics of hot, bulk QCD matter, which is expected to 
be created in ultra-relativistic heavy ion collisions at the SPS, the RHIC and 
the LHC. Our approach is capable of self-consistently calculating the freeze- 
out of the hadronic system, while accounting for the collective flow on the 
hadronization hypersurface generated by the QGP expansion. In particular, 
we perform a detailed analysis of the reaction dynamics, hadronic freeze-out, 
and transverse flow. 
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I. INTRODUCTION 



A major goal of colliding heavy-ions at relativistic energies is to heat up a tiny region 
of space-time to temperatures as high as are thought to have occured during the early 
evolution of the universe, a few microseconds after the big bang In ultra- relativistic 
heavy-ion collisions the four-volume of hot and dense matter, with temperatures above 
~ 150 MeV, is on the order of ~ (10 fm)^. The state of strongly interacting matter at such 
high temperatures (or density of quanta) is usually called quark-gluon plasma (QGP) p|. For 
a discussion of the properties and potential signatures of such a super-dense state see 

A particularly interesting aspect of producing such a hot and dense space-time region is 
that QCD, the fundamental theory of strong interactions, is expected to exhibit a transition 
to a new thermodynamical phase at a critical temperature Tq ~ 100 — 300 MeV. This 
phase transition has been observed in numerical studies of the thermodynamics of QCD at 
vanishing net baryon charge on lattices [0. It is the only phase transition of a fundamental 
theory that is accessible to experiments under controlled laboratory conditions. 

In this paper we shall investigate the dynamics of relativistic heavy ion collisions within 
a novel transport approach combining a macroscopic and a microscopic model. We shall 
focus here on collision systems currently under investigation at the CERN Super-Proton- 
Synchrotron (SPS), the Relativistic Heavy Ion Collider (RHIC) at BNL and the future Large 
Hadron Collider (LHC) at CERN. 

We shall work in natural units h = c = k = 1 throughout the paper. 

II. GENERAL ASPECTS OF MATCHING FLUID DYNAMICS TO 
MICROSCOPIC TRANSPORT 

In this section we discuss general aspects and assumptions of our model for the space-time 
evolution of high-energy heavy-ion reactions. In particular, we introduce fluid dynamics for 
the early, hot stage, and the matching to microscopic transport for the later, more dilute 
stages of the reaction. Within this section, quantities without subscript refer to the fluid, 
while properties of the microscopic transport theory carry the subscript micro. 
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A. Transport Equation for Incoherent Quanta / Particles 

The most basic assumption of our model for the evolution of high-energy heavy-ion 
reactions is that at the initial timeQ t = ti the highly excited space-time domain produced 
in the impact can be viewed as being populated by incoherent quanta on the mass-shell. 
Thus, the system can be described by a distribution function fi{x^,p'^), where x° = ti, 

= ^JfP + mf, and i labels different species of quanta. We will not discuss here how such a 
state of high entropy density could possibly be reached. That discussion is out of the scope 
of the present manuscript. Our work addresses the subsequent evolution of that initial state 
up to the so-called freeze-out of strong interactions in the system. 

The semi-classical evolution of the distribution function in the forward light-cone is 
described by means of a so-called transport equation, e.g. the Boltzmann equation |^ 

p-dMx^,pn = c, . (1) 

Ci is the collision kernel, describing gain or loss of quanta (particles) of species i in the phase- 
space cell around {x^^p") due to collisions. Note that we have dropped possible classical 
background fields in eq. (|I]). 

B. Moments of the Transport Equation: Hydrodynamics for the hottest stage 

For the problem at hand, however, the usefulness of eq. (|lD is rather limited. A major 
difficulty is that to obtain an analytical or numerical solution, in most cases one has to in- 
troduce an expansion of the collision kernel in terms of the number of incoming particles per 
"elementary" collision [^] (In most practical applications that expansion is even truncated 
at the level of binary collisions, 2 — > n). Obviously, the expansion is ill-defined at very 
high densities. A second major problem is to describe the hadronization process, i.e. the 
dynamical conversion of quarks and gluons into hadrons, on the microscopic level. Several 
interesting approaches to describe the hadronization of a plasma of quarks and gluons micro- 
scopically have been proposed in the literature, cf. e.g. Q and references therein. However, 
due to the very complicated nature of this process, many of those models have to involve 
some kind of ad-hoc prescriptions, which have quite significant impact on the results. A 
first-order QCD phase transition, as assumed in the following, is particularly difficult to 
model microscopically. 



Our choice of space-time variables is described in more detail below; for the moment, we assume 
that suitable variables have been chosen, and that the hypersurfaces of homogeneity are time- 
orthogonal everywhere. 
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At present we are not able to solve these problems in a fully satisfactory way. We can, 
however, circumvent them to some extent if we are mainly interested in the bulk dynamics 
of a hot QCD system. In this case we can employ relativistic ideal hydrodynamics [0] for 
the very dense stage of the reaction up to hadronization. 

Let us thus assume that it is feasible to employ the continuum limit. The first two 
moments of eq. (|ip yield the continuity equations for the conserved currents and for energy 
and momentum 

d-Ni = 0,d-e = 0. (2) 

In the following, we will explicitly consider only one conserved current, namely the (net) 
baryon current. All other currents, as e.g. strangeness, charm, electric charge etc., will be 
assumed to vanish identically (due to local charge neutrality and an ideal fluid) such that 
the corresponding continuity equations are trivially satisfied. 

Ideal fluid dynamics goes even further and assumes that the momentum-space distribu- 
tions in the local rest frame are given by either Fermi-Dirac or Bose-Einstein distribution 
functions, respectively. Dissipation and heat conduction which arise from higher moments 
are neglected. Since we restrict fluid-dynamics to the high-temperature and high-density 
stage, this approximation is at least logically consistent. In future work, it will be impor- 
tant though to check its quantitative accuracy. 

The density of secondary partons in the central region of high-energy nuclear collisions 
is very high. According to present knowledge, it is likely that the central region evolves 
from a stage of pre-equilibrium towards a QGP in local thermal equilibrium despite 
the large expansion rate. On the other hand, the very same calculations do not seem to 
support rapid chemical equilibration (in particular of the quarks), cf. also [Il2|-[T^. However, 
in most publications interactions among the secondary partons (and in particular particle 
production via inelastic processes) were treated perturbatively. Since the running coupling 
in a thermal plasma with T < 1 GeV is not very small, one can not exclude sizeable 



contributions from processes involving higher powers of |]I0[. Moreover, in addition to 
the semi-hard partons there might exist a coherent color field in the central region (between 
the receding nuclear "pancakes") which produces additional quark-antiquark pairs in its 



decay [15|. In any case, we will not argue in favor nor against rapid qq production and 
chemical equilibration but simply assume that the quark densities at the initial time of 
the hydrodynamical expansion are close to their chemical equilibrium values. At least for 
Pb+Pb at CERN-SPS energy, where experimental data already exists, this is basically the 
only way for our model to account for the fact that measured hadron multiplicity ratios are 
close to their chemical equilibrium values |]T6[. Since the expansion rate after hadronization 
is too large for "chemical cooking" (and in particular for strangeness equilibration), as 



will be discussed in section [IV D| , it would be virtually impossible to achieve approximate 
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chemical equilibrium during the later hadronic stages if starting from a QGP far off chemical 
equilibrium, cf. also [1T7| . 

The general picture as described above is summarized in the space-time diagram depicted 
in Fig. ^ We assume that ideal fluid dynamics is a reasonable approximation between 
the "initial" time ti and the hadronization hypersurface. After that, we will switch to a 
microscopic description employing the binary collision approximation for the collision kernel. 
In particular, we will employ the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) 
transport model, see below. 



C. Microscopic Transport from Hadronization to Freeze-Out 



One may ask why it is not sufficient to rely on hydrodynamics up to some rather late 
stage of the reaction, after which one postulates that all particle momenta are "frozen" and 
thus are equal to those measured in the detector at t/ — * oo. That approach has been 
applied to nuclear collisions by many authors, for recent work cf. e.g. |T^^I|, and leads to 
reasonable results for the single-particle spectra of the most abundant hadron species vr, K, 
p, A. However, the following limitations arise: 

as 



First, the evolution must clearly be non-ideal in the late stages of the reaction |22 



the system approaches "freeze-out" . This can manifest in decoupling of various components 
of the fluid (e.g. pions and nucleons), i.e. each component develops an individual collective 
velocity [^. Another aspect is that the —pdV expansion work performed by the fluid can 
be partly compensated by entropy production (+TdS') such that the expansion may even 
become isoergic, dE = 0, instead of isentropic, dS' = |^ . 

Moreover, since each hadron is propagated individually, and its interactions with other 
hadrons are described on the basis of elementary processes, microscopic transport models 
offer the opportunity to calculate the freeze-out conditions instead of just putting them in by 



hand as is done in the purely fluid-dynamical approaches p8|-pi[]. There one assumes that 
freeze-out occurs whenever some criterion is fulfilled, e.g. when the temperature drops below 
some "guessed" value. In contrast, the non-truncated transport eq. (|T|) can describe self- 
consistently the freeze-out of the system: no decoupling hypersurface is imposed by hand, 
but rather is determined by an interplay between the (local) expansion scalar d ■ u [||,^^-|27 



(where u is the four- velocity of the local rest-frame), the relevant elementary cross sections 
and decay rates, and the equation of state (EoS), which actually changes dynamically as 
more and more hadron species decouple. This is obviously a key point for being able to study 
and predict the dependence of the final state on collision energy (i.e. on the initial entropy 
or energy density), system size etc., instead of just fitting it by an appropriate choice of a 
freeze-out hypersurface. Note e.g. that the nucleons emerging from the QCD hadronization 
phase transition in the early universe were able to maintain chemical equilibrium down 
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to temperatures of about ~ 50 MeV In heavy-ion collisions at CERN-SPS energies, 
however, one finds chemical freeze-out temperatures on the order of 140 — 160 MeV ||16|| . 
The origin of this difference lies in the much smaller expansion rate (Hubble constant) of 
the early universe as compared to a high-energy heavy ion collision |^^, and can only be 
explained within kinetic theory but not within pure hydrodynamics. 

Another complication arises from the fact that close to the freeze-out hypersurface the 



freeze-out process feeds back on the evolution of the fluid [28-30]. This will in general 
deform the freeze-out hypersurface, say an isotherm of given temperature Tfo- It will differ 
from that found a posteriori from the solution of eqs. (Q) in the whole forward light-cone. 
Furthermore, the idealization that the transition from ideal flow to free-streaming occurs on 
a sharp hypersurface, i.e. a three- volume in space-time, is rather crude. One instead expects 
a smooth transition as the temperature (and the density of particles) decreases, cf. e.g. the 



discussion in [31]. This is supported by studies of the hadron kinetics close to freeze-out 



with realistic cross-sections ]|32| , |33| , cf. also section [IV C . 

Finally, it is likely that the freeze-out is not universal for all hadron species, simply 
because their transport cross-sections are very different. One can therefore hardly assume 
that all hadron species decouple on the same hypersurface [I26|j3^ -p4| . The clearest exam- 
ple for this is the transverse momentum distribution of Q baryons obtained by the WA97 



collaboration [^ for Pb+Pb collisions at CERN-SPS energy, ^/s = 17 A GeV. Unlike is the 



case for pions, nucleons, anti-nucleons, and lambdas, the pr distribution of omegas as calcu- 



lated within hydrodynamics with freeze-out on the T = Tfo = 130 MeV hypersurface [^] is 
much stiffer than the experimental flnding. Indeed, more detailed kinetic treatments which 
explicitly account for the small transport cross-section of Q baryons in a meson-rich hadron 
gas, emerging either from fragmentation of longitudinally stretched color-strings [|36| or an 



incoherent hot plasma of quarks and gluons [37], show that these multiple-strange particles 



freeze out earlier and pick up less collective transverse flow than pions and nucleons, for 
example. 

D. Transition from Fluid Dynamics to Microscopic Transport 

A few remarks on the transition from hydrodynamics to microscopic transport are in 
order here. In general, one should introduce source terms —d ■ Nmicro and —d ■ Qmicro on the 
right-hand-sides of eqs. (^, where 

^microi-^) ~ J ^^/i, micro (3^; ^) ; (3) 
^ microi.-^) ~ j jjy k^k /j, micro (2^; k) (4) 
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denote the net baryon current and the energy-momentum tensor of the microscopic transport 
model, respectively. Accordingly, external sources of particles have to be introduced in the 
transport equation, which model the net-baryon charge and energy-momentum transfer from 
the fluid. This way a self-consistent solution in the whole forward light-cone, starting from 
the initial hypersurface t = U could be obtained. 

However, if a space-time region bounded by a hypersurface cx^ exists where fluid- 
dynamics is an adequate approximation, one can choose an arbitrary hypersurface cr'^wUch 
within this region where to switch from eqs. (|]) to (|l]). One can then simply assume 
N^icro = and Qmicro = in the interior^], and = 0, = in the exterior. On 
that hypersurface, one sets Nmicro = N, Qmicro = ©• This is because hydrodynamics is a 
limiting case of eq. (|lD, and this more general transport equation will automatically recover 
the fluid-dynamical solution in the space-time region between (Jg^nch ^h- 

For the particular model discussed here, cr^^^^^^ can not precede the hadronization hy- 
persurface since our microscopic transport model deals with color-singlet states, only. Also, 
it employs the binary collision approximation of the kernel, which becomes less justified in 
the hot and dense stage preceding hadronization. 

Furthermore, as will be discussed in more detail below, in high-energy heavy-ion collisions 
it turns out that the boundary of validity of (ideal) fluid dynamics, cr^, can not extend far 
into the post-hadronization stage (the less the higher the collision energy). Thus, we conclude 
that the hadronization hypersurface is the most natural choice for the switch from eqs. @) 
to (0). 

The phase-space distribution of particles of species i on cr'^y^^^^^ is then given by |^ 

Ei^^ = Jda-p f{p-u), (5) 

where u'^ is the four- velocity of the local rest-frame. An explicit expression for the geometry 
suitable for high-energy collisions will be given below. For time-orthogonal hypersurfaces as 
depicted in Fig. |l| one has 



daf. 



= {d'x, 0) , (6) 

t=const. 



/ i,micro' 



the left-hand-side of eq. (^ being simply Eifi^ 

It is clear that by construction the microscopic transport starts from a state of local 
equilibrium on (j'^yjuchi ^^e hypersurface where the switch is performed. The local energy 
density, net baryon density, and collective expansion velocity are those obtained from the 
hydrodynamical solution. Thus, the conserved currents and the energy-momentum tensor 
of the microscopic transport theory assume the form appropriate for ideal fluids 0, 



2 a- 



Interior" meaning towards the origin of our space-time diagram, Fig. |l|. 
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^micro , (7) 

^micro ~^ Pmicro) U^U Pmicrod^ • (^) 

Now, in order that Qmicro = on (Jg^nchy pressure at given energy and baryon density 
must equal that of the fluid-dynamical model, i.e., the equations of state in local thermody- 
namical equilibrium must be the same. In general this requirement is non-trivial. For ideal 
gases, however, it can be obeyed by simply including the same states i in the microscopic 
transport (0) as in the grand partition function which is used to calculate the equation of 
state employed in hydrodynamics. We shall discuss this point in more detail when presenting 
our specific equation of state below. 

We finally briefly discuss one last aspect of the switch from fluid dynamics to microscopic 
transport on some hypersurface cr switch- -^^ already mentioned above, this hypersurface is 
assumed to be within the region of validity of ideal hydrodynamics, and should be identified 
with the hadronization hypersurface. However, the latter will in general also exhibit time- 
like parts (points where the normal vector on Cg^nch space- like). A schematic example 
is given in Fig. ^ The initial condition of the microscopic transport on o"^^^^^^ can now 
not be chosen arbitrarily. This is clear from the fact that the points 1 and 2, for example, 
are causally connected. The simplest way to prevent violation of the evolution equations is 
to specify initial conditions on a purely space-like hypersurface (e.g. t = ti) and to employ 
the dynamical equations, in our case the continuity equations (||), to calculate 6 and on 
^switch- This way the states of the system at 1 and at 2 are consistent. 

However, one problem with switching on a hypersurface with time-like parts remains. 
As discussed above, eq. (^) conserves energy-momentum and (net) baryon charge. For this 
to hold, it actually does not count the flow of the currents from the inside to the outside of 
^'switch but, actually, the net flow. That is, the difference of outflowing and inflowing charge, 
momentum etc. The inflow is due to those parts of the thermal distribution function f{x,p) 
which move into the opposite direction than the fluid. Due to the exponential tails of / such 
particles clearly always exist, but their number decreases strongly if the collective flow is 
strong. In this case the locally isotropic momentum-space distribution is strongly boosted. 

Thus, the in-current is obtained under the assumption that within an infinitesimal region 
on both sides of the hypersurface there is hydrodynamic flow and local thermodynamical 
equihbrium. For this reason, c'^^^^^f^ must be entirely within the region of validity of fluid 
dynamics, cr^. Again, in this case eq. (^ gives the net flow of all the currents from the 
fluid-region to the region where we apply the microscopic transport. 

The problem is, however, that in some part of momentum and coordinate space the left- 
hand-side of eq. (^ can be negative. This means that the ingoing flow exceeds the outgoing 



flow. These "negative contributions" were already discussed by several authors [28,29,39 



Since we will interpret d^pda ■ pf{p ■ u)/E as a probability distribution we have to require 
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positive definiteness. This can either be achieved by multiplying with a cut-off function 
0(d(T ■ p), which leads to a slight violation of the conservation laws; or by integration over 
sufficiently large bins in momentum and coordinate-space, and random redistribution of the 
particles within the bins, which smears out the distribution over momentum and coordinate 
space. 

A rigorous solution of this problem requires to introduce the above-mentioned source- 
terms in the fluid-dynamical evolution equations as well as in the microscopic transport 
equation. However, for the cases studied here the negative contributions were not relevant. 
The main reason is that the collective flow velocity on the time-like parts of the hyper- 
surface is close to one, such that net flow of particles from the microscopic transport to 
hydrodynamics does not occur. For non-relativistic flow on 0"^^^^^^, however, the negative 
contributions would be more serious. 



III. SPECIFIC MODEL FOR HIGH-ENERGY HEAVY-ION COLLISIONS 



In the present paper we shall use hydrodynamics to model a first order phase transition 
from a QGP to a hadronic fluid, and combine it with a microscopic transport calculation 
for the later, purely hadronic stages of the reaction. 

In the following sections we describe the particular hydrodynamical and transport models 
employed here, cf. also PTI and PD|. 



A. Scaling Hydrodynamics 

As already mentioned above, hydrodynamics for hadronic collisions is defined by (local) 
energy-momentum and net baryon charge conservation, 

d-e = 0,d-NB = 0. (9) 

Q'^" denotes the energy-momentum tensor, and A'^^ the current of net baryon charge. 

For ideal fluids, the energy-momentum tensor and the net baryon current assume the 
simple form |^ 

= {e + p) u^'u" - pg"" , A^^ = peu^" , (10) 

where e, p, pB are energy density, pressure, and net baryon density in the local rest frame 
of the fluid, which is defined by A^^ = (pB,C)). Let us, in the following, work in the metric 
g^'' = diag(+, — , — , — ). = 7(1, v) is the four-velocity of the fluid {v is the three-velocity 
and 7 = (1 — t?^)~^/^ the Lorentz factor). The system of partial differential equations (P) is 
closed by choosing an equation of state (EoS) in the form p = p{e,pB), cf. below. 
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For simplicity, we assume a cylindrically symmetric transverse expansion with a longitu- 
dinal scaling flow profile, Vz = z/t At z = 0, equations (||) reduce to 

dtE + dT[{E+p)vT] =-(^ + i)(^+J9) , 

dtM + driMvT+p) =-(— + ]) M , (11) 

dtR + Ot (Rvt) = - + ^) i? , 

where we defined E = 0°°, M = Q^^, and R = N^. In the above expressions, the index T 
refers to the transverse component of the corresponding quantity. 

The set of equations ([TT|) describes the evolution in the z = plane. Due to the assump- 
tion of longitudinal scaling, the solution at any other z ^ can be simply obtained by a 
Lorentz boost. The above equations also imply 



dp 
dr] 



, (12) 



T,rT 



where r] = Artanh Vz and r = \fW—z^. This means that on r = const, hypersurfaces 
pressure gradients in rapidity direction vanish, and there is no flow between adjacent in- 
finitesimal rapidity slices. However, only for net baryon free matter, pB = 0, does this 
automatically also mean that the temperature T is independent of the longitudinal fluid 
rapidity r]. In the case 7^ equation (|1^) only demands 



dT 
dr] 



+ pB 



dr] 



. (13) 



s and ]iB denote entropy density and baryon-chemical potential, respectively. If other charges 
like strangeness or electric charge are locally non-vanishing, additional terms appear. Equa- 
tion (|13D does not imply that the rapidity distribution of produced particles is flat (i.e. 
independent of rapidity) or that the rapidity distributions of various species of hadrons, 
e.g. pions, kaons, and nucleons, are similar. Any rapidity- dependent T and /i^ that satisfy 
eq. (0) are in agreement with energy-momentum and net baryon number conservation, 
as well as with longitudinal scaling flow Vz = z/t Note also that non-trivial solu- 

tions of eq. (0) in general also yield dfis/dr] 7^ on the hadronization hypersurface (i.e. a 
rapidity-dependent strangeness-chemical potential), even if the strangeness-density ps = 
everywhere in the forward light-cone. In this paper, however, we do not explore the rapidity 
dependence of the particle spectra, and thus simply assume that T and pb are independent 
of 77. 

The fluid-dynamical evolution equations can be solved numerically on a discretized space- 
time grid, cf. e.g. [^ , ^ . 
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B. Equation of State 



To close tlie system of coupled equations of liydrodynamics, an equation of state (EoS) 
lias to be specified. From eq. (|l^) it follows that for an ideal gas the pressure p is given by 

p = q.{q.Q) , (14) 

where is orthogonal to u'^ and normalized to g ■ g = — 1. In particular, in the local rest- 
frame uf" = (1, 0), we can choose qf" oc (0, 1, 0, 0) + (0, 0, 1, 0) + (0, 0, 0, 1). Then, from the 
definition of the energy- momentum tensor from kinetic theory, eq. (^), we obtain 

p(T,/iB,/is) = ^y"^y/i(fc;T,/ii3,/i5) . (15) 

The sum over i extends over the various particle species. The grand canonical potential is 
given hj Q = —pV, where V = J da ■ u denotes the three-volume of the given hypersur- 
face of homogeneity. All other quantities can be obtained via standard thermodynamical 
relationships. E.g., the densities of entropy, net baryon charge, and energy are given by 

s{T,fiB,f^s) = , (16) 

Pb{T,i2b,I^s) = ^ , (17) 

duB 

Ps{T,ijb,I^s) = ^ = , (18) 

e{T,pB,fJ's) = Ts-p + pbPb ■ (19) 

From p(T, pb, f^s), Pb(T, pb, fJ's) and e(T, pb, fJ's) one can construct the function p{e, pb) 
which is needed to close the system of continuity equations (H). 

So far, we discussed an ideal gas, only. However, lattice QCD predicts a phase transi- 
tion from ordinary nuclear matter to a so-called quark-gluon plasma (QGP) at a critical 
temperature of Tc = 140 — 160 MeV (for pB = 0). We will employ a very simple and 
intuitive, though not very well justified description of this phase transition. We model the 
high-temperature phase as an ideal gas of u, d, s quarks (with masses tjIu = ma = 0, 
TUs = 150 MeV), and gluons, employing the well-known MIT bag model EoS P,^. In this 



model the non-perturbative interactions of the "deconfined bag" of quarks and gluons with 
the true vacuum are parameterized by a bag constant B. To make this state thermody- 
namically unfavorable at low temperatures, the bag contribution to the pressure must be 
negative. Thus, when computing the pressure of the QGP phase we subtract B from the 
right-hand-side of eq. (^). Accordingly, the energy density receives a positive contribution, 
cf. (0), while s and pb remain unchanged. This additional "bag term" can also be under- 
stood as an additional contribution -j-Bg'^" due to the non-perturbative interactions to the 
energy-momentum tensor of the QGP-fiuid. 
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In the low-temperature region we assume an ideal hadron gas that includes the well- 
established (strange and non-strange) hadrons up to masses of ~ 2 GeV. They are listed in 
Tab. I and Although heavy states are rare in thermodynamical equilibrium, they have a 
larger entropy per particle than light states, and therefore have considerable impact on the 
evolution. In particular, hadronization is significantly faster as compared to the case where 
the hadron gas consists of light mesons only (see the discussion in |r8|P^ , pT] , |50| , ^ ) . 

The actual model used for the hadronic stage of the reaction (UrQMD, see section 111 D ) 
additionally assumes a continuum of color-singlet states called "strings" above the m ~ 
2 GeV threshold to model 2 ^ n processes and inelastic processes at high CM-energy. For 
example, the annihilation of an p on an Q is described as excitation of two strings with 
the same quantum numbers as the incoming hadrons, respectively, which are subsequently 
mapped on known hadronic states according to a fragmentation scheme. Since we shall be 
interested in the dynamics of the fi-baryons emerging from the hadronization of the QGP, 
it is unavoidable to treat string-formation. The fact that string degrees of freedom are not 
taken into account in the EoS (|T^) does not represent a problem in our case because we 
focus on rapidly expanding systems where those degrees of freedom can not equilibrate E? 



The phase coexistence region is constructed employing Gibbs' conditions of phase equi- 
librium. The bag parameter of -B = 380 MeV/fm^ is chosen to yield the critical temperature 
Tc ~ 160 MeV at pb = 0. By construction the EoS exhibits a first order phase transition, 
as is also expected in QCD for the quark-hadron phase transition in the case of three light 
quark flavors |^ . 

The most striking aspect of a first-order phase transition with respect to the dynamical 
evolution is that the pressure is almost constant within the phase coexistence region (in fact, 
in a fluid where all conserved currents vanish identically p = const, within the mixed phase). 
Thus, the isentropic speed of sound. 



2 dp 



, , (20) 

s/pB 



is very small. This quantity characterizes the pressure gradient caused by a given energy 
density gradient along an isentrope, i.e. at constant entropy density per net baryon density 
(recall that all continuous solutions of the relativistic ideal-fluid dynamical equations con- 
serve the entropy). A very small Cs means that (isentropic) expansion is inhibited because 
the fluid does not "respond" to energy density gradients. In heavy-ion collisions this reflects 
in a particularly "soft" expansion if the mixed phase occupies the largest space-time volume 
of all three phases [pl] , ^ ,^g|-|5T| . For recent discussions of the consequences of this effect in 
cosmology (primordial black hole formation, evolution of density perturbations through the 
QCD phase transition) see e.g. . 

However, a nearly vanishing isentropic velocity of sound does only occur if the net baryon 
density is not very large, as e.g. in the cosmological QCD phase transition or in the central 
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region of high-energy colhsions studied here. In heavy-ion colhsions at much lower energies, 
where the net baryon density in the central region is rather large, cs is not very small. 
Despite the first-order phase transition, the isentropic expansion of baryon-dense fluids is 
not inhibited |[53|] . 



One should also be aware of the fact that by constructing the phase coexistence region 
with Gibbs' conditions we implicitly assume a "well-mixed" phase, i.e. that the transition 
from the QGP to the hadronic stage proceeds in equilibrium. This is the common approach 



widely employed in the literature [||^,||,||,||,|3,||,||,||,|9|^j5|^ and so far it is 
not in contradiction to existing data. It is based on the picture that the first-order phase 
transition proceeds via nucleation of hadronic bubbles in the expanding QGP [Q, and that 



the bubble nucleation and growth is fast as compared to the expansion rate such that the 
two phases are approximately in pressure equilibrium. However, this scenario is less likely 
to apply to high-energy heavy-ion collisions than to the cosmological QCD phase transition. 



because in the former case the expansion rate is many order of magnitude larger [27|. In 



particular, it has been speculated recently that the time-scale for supercooling down to the 



spinodal instability is comparable to that for homogeneous bubble nucleation ||5j. Thus, 
it may well be that the phase transition proceeds via spinodal decomposition rather than 
bubble nucleation. In that case, the "soft" mixed phase with c^^ ^ Q would be absent and 
shorter reaction times may be expected. In any case, we postpone a detailed dynamical 
study of the latter scenario to a future publication, and shall restrict ourselves here to the 
more conservative picture assuming an adiabatic phase transition. 
Finally, we have to specify the initial conditions: 

SPS: For collisions at SPS energy we assume that hydrodynamic flow sets in on the hyperbola 
Tj = 1 fm/c. This is a value conventionally assumed in the literature, cf. e.g. PT| . 



We further employ a (net) baryon rapidity density (at mid-rapidity) of dNB/dy = 
80, as obtained by the NA49-collaboration for central Pb-|-Pb reactions The 
average specific entropy in these collisions is s/p^ = 45±5 (the bar indicates averaging 
over the transverse plane). That entropy per net baryon fits most measured hadron 



multiplicity ratios within ±20% |jT6|. The corresponding initial energy and net baryon 
densities (e, = 6.1 GeV/fm^, = 4.5po) are assumed to be distributed in the transverse 
plane according to a so-called "wounded nucleon" distribution with transverse radius 



Rt = 6 fm, i.e. e{Ti,rT), PB{Ti,rT) oc /(rr), with /(ry) = fyl - r'^/Rj-- The initial 
temperature and quark-chemical potentials are (they are of course not exactly constant 
over the transverse plane) Tj ~ 220 MeV, Pq ^ 150 MeV, Ps = 0. The transverse 
velocity field on the t = Ti hyperbola is assumed to vanish. 

RHIC: Due to the higher parton density at mid-rapidity as compared to collisions at SPS 
energy, thermalization may be reached earlier at RHIC. According to various stud- 
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ies [§,0, thermalization might occur within ~ 0.3 — 1 fm. We assume Xj = Rt/IQ = 
0.6 fm. The net baryon rapidity density and specific entropy at mid-rapidity in central 
Au+Au at ^/s = 200A GeV is predicted by various models of the initial evolution, e.g. 
the parton cascade model PCM, RQMD 1.07, FRITIOF 7, and HIJING/B, to be in the 
range dNs/dy ^ 20 - 35, s/pb ^ 150 - 250 0,13 . We will employ dNs/dy = 25 
and s/pb = 205 li = 20 GeV/fm'^, p^ = 2.3po)- These parameters could of course 
be fine-tuned once the first experimental data are available. As in the above case, 
e(Tj,TT) and PBijij^T) are initially distributed in the transverse plane according to a 
wounded-nucleon distribution with Rt = 6 fm. The initial temperature and quark- 
chemical potentials follow as Tj ~ 300 MeV, ~ 45 MeV, fig = 0, respectively. This 
corresponds to a transverse energy on the t = Ti hyperbola of dEx/dy ^1.3 TeV, 
which decreases to dEx/dy ~ 720 GeV on the hadronization hypersurface . 



LHC: The initial conditions for CERN-LHC energy are, of course, less well known. Qualita- 
tively, and according to present expectations, it appears reasonable to assume that 

1. the density of minijets produced at time tq ~ l/po, where po ~ 2 GeV is the 
minijet cut-off scale, is much larger than at BNL-RHIC energy. The most recent 
estimates of the energy densities in the central region span the range eo = (0.3 — 
1.3) TeV/fm^ |]6T|j6^ . The results to be expected from p + p + A, and A + A 



at BNL-RHIC will probably not reduce the uncertainties by much because the 
energy density at y ~ and y/s = 5.5A TeV depends strongly on the model for 
the nuclear parton distribution functions at very small x, out of range for RHIC. 

2. the higher initial density of partons could also lead to somewhat faster equili- 
bration than at the lower energies. Note that the produced gluons already have 
the "right" thermal energy per particle, eo/po — 2.7To [0. The distribution in 
momentum-space, however, has to become isotropic via rescattering among the 
partons [0,0. 

3. the net baryon charge in a rapidity-slice Ay = 1 around y = is even smaller than 
at RHIC, and can in practice be neglected if one focuses on the bulk dynamics of 
the central region (in the same way as we neglect net strangeness, charm, etc.). 
Note however, that smaller rapidity bins may exhibit quite large fluctuations of 
the initial net baryon charge |]63|. 



Thus, in view of these uncertainties, it is clear that precise quantitative predictions 
for the CERN-LHC energy are hardly possible at the moment. Our more modest aim 
will therefore be to discuss a set of even more "extreme" initial conditions than those 
employed for CERN-SPS and BNL-RHIC energies, to give an idea how the dynamical 
evolution may continue at even higher energies. Whether or not that set of initial 
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conditions corresponds closely to the LHC case can not be decided presently on solid 
grounds. 

Thus, we employ a thermalization time = 0.3 fm, an initial energy density (not 
on the To but on the Xj hypersurface !) q = 230 GeV/fm^, and a vanishing net 
baryon charge, dA^'^/d?/ = pi = 0. Again, the initial energy density is distributed in 
the transverse plane according to a wounded nucleon distribution with Rt = 6 fm. 
The initial temperature is about Tj ~ 580 MeV (it is not exactly constant over the 
transverse plane), the initial transverse energy is dEx/dy ~ 7.8 TeV. 



C. Hadronization and the transition to microscopic dynamics 



Having specified the initial conditions on the r = Tj hypersurface and the EoS, the hydro- 
dynamical solution in the forward light-cone is determined uniquely. As already mentioned 
in section P D| , we assume that it is not a bad approximation to determine the hadroniza- 
tion hypersurface a posteriori from the solution in the whole forward light-cone. In other 
words, the hadronization hypersurface is assumed to be within the region of validity of 
hydrodynamics. 



In parametric representation, the hypersurface is a function of three parameters 



In our case, due to the symmetry under rotations around and Lorentz-boosts along the beam 
axis, two of these parameters can simply be identified with rj and 0, while r and tt depend 
only on the third parameter, call it (. Thus, ( G [0, 1] parameterizes the hypersurface in the 
planes of fixed rj and (in the mathematically positive orientation, i.e. counter clock- wise). 
The normal is 0] 

9cr" 9(7^ da'^ _ , , , 
da, = W7^^^dCdr;d0 

/ dvT , dr dr . dr^ . , \ i /■ i i , /r.-i \ 

cosh 7], — cos<p, — sm0, —— smhrj ttt d(^drid(p . [21) 



V dC " dC dC dC 

This expression naturally looks simpler in the (r, r], r-p, 0)-basis, cf. e.g. but we will 
nevertheless write all vectors and tensors in the (t, x, y, z)-basis throughout the manuscript, 
even if the components are written in terms of the variables r, rj, r^, and 0. 

We can now apply eq. (^ to compute the number of hadrons of species i hadronizing at 
space-time rapidity rj, proper time r, and position rx (cos (x ~ 0) ? sin (x — 0)), with four- 
momentum = {rriT cosh y, px cos x, Pt sin x, "^t sinh y) , 

ttt \ pt cos (x - 0) -77 - cosh(?/ - rj) —— ] Ji[p- u) . (22) 



d'^prdydrfd^dcf) \ d( d( 

= 7r(cosh 7], vt cos(x — 4>),vt sin(x — 0), sinh rj) denotes the fiuid four- velocity. Thus, the 
direction of the particle momentum in the transverse plane is determined by the angle X; 
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while the relative angle between and the transverse flow velocity, vt-, is denoted by 0. 
/ is either a Bose-Einstein or Fermi-Dirac distribution function, depending on the particle 
species under consideration. 

How is the distribution (^) actually passed to the microscopic model ? First, it is 
integrated over space-time (17, 0) and momentum space (pr, y), rounded to an integer 
value (the hadronic transport model described in the next section deals with integer number 
of particles, only), and the distribution (^) divided by Ni is used as probability distribution 
to randomly generate space-time and momentum-space coordinates for Ni hadrons of species 
i. Of course, due to the fact that our system has a surface and does not extend to inflnity in 
the transverse plane, hadronization does not occur on a r = const, hypersurface, cf. Fig. |^. 
Thus, if we look at our expanding system on r = const, surfaces, there exists an interval 
where the two models, hydrodynamics and the microscopic transport, are applied in parallel. 



D. Microscopic dynamics: the UrQMD approach 

The ensemble of hadrons generated accordingly is then used as initial condition for 
the microscopic transport model Ultra-relativistic Quantum Molecular Dynamics (UrQMD) 
J^. The UrQMD approach is closely related to hadronic cascade |6^, Vlasov-Uehling- 



Uhlenbeck [^| and (R)QMD transport models |^|. We shall describe here only the part 



of the model that is important for the application at hand, namely the evolution of an 
expanding hadron gas in local equilibrium at a temperature of about Tc ~ 160 MeV. 
The treatment of high-energy hadron-hadron scatterings, as it occurs in the initial stage of 
ultrarelativistic coUisions, is not discussed here. A complete description of the model and 
detailed comparisons to experimental data can be found in pOl . 

The basic degrees of freedom are hadrons modeled as Gaussian wave-packets, and strings, 
which are used to model the fragmentation of high-mass hadronic states via the Lund scheme 
p8| . The system evolves as a sequence of binary collisions or 2 — A^-body decays of mesons, 
baryons, and strings. 

The real part of the nucleon optical potential, i.e. a mean-fleld, can in principle be 
included in UrQMD for the dynamics of baryons (using a Skyrme-type interaction with a 
hard equation of state). However, currently no mean fleld for mesons (the most abundant 
hadrons in our investigation) are implemented. Therefore, we have not accounted for mean- 
fields in the equation of motion of the hadrons. To remain consistent, mean fields were 
also not taken into account in the EoS on the fiuid-dynamical side. Otherwise, pressure 
equality (at given energy and baryon density) would be destroyed. We do not expect large 
modifications of the results presented here due to the effects of mean fields, since the "fiuid" 
is not very dense after hadronization and current experiments at SIS and AGS only point 
to strong medium-dependent properties of mesons (kaons in particular) for relatively low 
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incident beam energies {Eiab < 4 GeV/nucleon) [^. Nevertheless, mean fields will have to 
be included in the future; a fully covariant treatment of baryon and meson dynamics within 
UrQMD derived from a chiral Lagrangian []7D[ is currently under development. 

Binary collisions are performed in a point-particle sense: Two particles collide if their 
minimum distance d, i.e. the minimum relative distance of the centroids of the Gaussians 
during their motion, in their CM frame fulfills the requirement: 



d <do = Y o-tot = o-{y/s, type). (23) 

The cross section is assumed to be the free cross section of the regarded collision type {N—N, 
N - A, TT - N ...). 

The UrQMD collision term contains 53 different baryon species (including nucleon, delta 
and hyperon resonances with masses up to 2 GeV) and 24 different meson species (including 
strange meson resonances), which are supplemented by their corresponding anti-particle and 
all isospin-projected states. The baryons and baryon-resonances which can be populated in 
UrQMD are listed in table |, the respective mesons in table y - full baryon/antibaryon 
symmetry is included (not shown in the table), both, with respect to the included hadronic 
states, as well as with respect to the reaction cross sections. All hadronic states can be 
produced in string decays, s-channel collisions or resonance decays. 

Tabulated and parameterized experimental cross sections are used when available. Res- 
onance absorption, decays and scattering are handled via the principle of detailed balance. 
If no experimental information is available, the cross section is either calculated via an One- 
Boson-Exchange (OBE) model or via a modified additive quark model which takes basic 
phase space properties into account. 

In the baryon-baryon sector, the total and elastic proton-proton and proton-neutron 
cross sections are well known 0] . Since their functional dependence on ^/s shows a compli- 
cated shape at low energies, UrQMD uses a table-lookup for those cross sections. However, 
many cross sections involving strange baryons and/or resonances are not well known or even 
experimentally accessible - for these cross sections the additive quark model is widely used. 

As we shall see later, the most important reaction channels in our investigation are 
meson-meson and meson-baryon elastic scattering and resonance formation. For example, 
the total meson-baryon cross section for non-strange particles is given by 

(^M^iVs) = J2 Ob' mBjM, mnW Jr, Mr) ^'^vtc^ . r\ 

R=A,N* l^'^'B + -'-Jl^'^M + ij 

^R^MB^tot 

X ^ ^ (24) 

PCMS (Mfi - v^)2 + %i 

with the total and partial -y/i-dependent decay widths Ttot and Th^mb- The full decay 
width Ttot{M) of a resonance is defined as the sum of all partial decay widths and depends 
on the mass of the excited resonance: 



17 



TtotiM) = J2 r^,,(M) . (25) 

br={i,j} 

The partial decay widths Tij{M) for the decay into the final state with particles i and j is 
given by 

here Mr denotes the pole mass of the resonance, T^^ its partial decay width into the channel 
i and j at the pole and / the decay angular momentum of the final state. All pole masses 



and partial decay widths at the pole are taken from the Review of Particle Properties |7I 
TijlM) is constructed in such a way that TijlMn) = F^-' is fulfilled at the pole. In many 
cases only crude estimates for F^-' are given in [|7T| - the partial decay widths must then 
be fixed by studying exclusive particle production in elementary proton-proton and pion- 
proton reactions. Therefore, e.g., the total pion-nucleon cross section depends on the pole 
masses, widths and branching ratios of all A^* and A* resonances listed in table |I[ Resonant 
meson-meson scattering (e.g. tt + ti^potti + K-^ K*) is treated in the same formalism. 

In order to correctly treat equilibrated matter (we repeat that the hadronic matter 
with which UrQMD is being initialized in our approach is in local chemical and thermal 
equilibrium), the principle of detailed balance is of great importance. Detailed balance is 
based on time-reversal invariance of the matrix element of the reaction. It is most commonly 
found in textbooks in the form: 

^ Pi 9i ^ /Qr,\ 

^/-» = ' ^27) 

Pf 9f 

with g denoting the spin-isospin degeneracy factors. UrQMD applies the general principle 
of detailed balance to the following two process classes: 

1. Resonant meson-meson and meson-baryon interactions: Each resonance created via a 
meson-baryon or a meson-meson annihilation may again decay into the two hadron 
species which originally formed it. This symmetry is only violated in the case of three- 
or four-body decays and string fragmentations, since N-body collisions with (N> 2) 
are not implemented in UrQMD. 

2. Resonance- nucleon or resonance-resonance interactions: the excitation of baryon- 
resonances in UrQMD is handled via parameterized cross sections which have been 
fitted to data. The reverse reactions usually have not been measured - here the princi- 
ple of detailed balance is applied. Inelastic baryon-resonance de-excitation is the only 
method in UrQMD to absorb mesons (which are bound in the resonance). Therefore 
the application of the detailed balance principle is of crucial importance for heavy 
nucleus-nucleus collisions. 
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Equation (pTl), however, is only valid in the case of stable particles with well-defined 
masses. Since in UrQMD detailed balance is applied to reactions involving resonances with 
finite lifetimes and broad mass distributions, equation ( P7D has to be modified accordingly. 
For the case of one incoming resonance the respective modified detailed balance relation has 



been derived in Here, we generalize this expression for up to two resonances in both, 
the incoming and the outgoing channels. 

The differential cross section for the reaction (1 , 2) (3 , 4) is given by: 

here the pi in the 5-function denote four-momenta. The 5-function ensures that the par- 
ticles are on mass-shell, i.e. their masses are well-defined. If the particle, however, has a 
broad mass distribution, then the 5-function must be substituted by the respective mass 
distribution (including an integration over the mass): 

da?^ = - — ^ dn Y\pu 9 ■ (29) 

Incorporating these modifications into equation (|^) and neglecting a possible mass- 
dependence of the matrix element we obtain: 



d4l _ {pi,) (2g, + l)(2g, + l) ^ dalt 

(2^3+1) (2^4+1) h. ^^^"^^^^"^^11-^^^^ ^ ■ 



Here, Si indicates the spin of particle i and the summation of the Clebsch-Gordan-coefficients 
is over the isospin of the outgoing channel only. For the incoming channel, isospin is treated 
explicitly. The summation limits are given by: 

J_ = max ( I ji - j2 1 , lis - J4 1 ) (31) 
J+ = min (ji + ja, js + ji) ■ (32) 

The integration over the mass distributions of the resonances in equation (^U|) has been 
denoted by the brackets (), e.g. 



4 , 

with the mass distribution A^{m) given by a free Breit-Wigner distribution with a mass- 
dependent width according to equation (pS]): 

1 r(m) 

AJm) = —-. -TT — — — rrr-— with \im Ar(m) = 5 (rrir — m) , (33) 

with the normalization constant 
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N = I 7 ^f^^L ,,,, dm . (34) 

' {rur - m)2 + r(m)V4 ^ ^ 



-oo 



Alternatively one can also choose a Breit-Wigner distribution with a fixed width, the nor- 
malization constant then has the value = 27i. 

The most frequent applications of equation ( pOD in UrQMD are the processes A1232 
A^A^ and A1232 A1232 ^ A^A^. 



IV. RESULTS FOR HEAVY-ION COLLISIONS AT CERN-SPS, BNL-RHIC AND 

CERN-LHC 



We now present some representative results for central collisions of heavy ions at CERN- 
SPS, BNL-RHIC and CERN-LHC energies. We will focus on single-inclusive momentum- 
space distributions and the space-time picture of freeze-out following the hadronization phase 
transition. As we shall see, one already gains much insight into the dynamics of high-energy 
heavy-ion collisions from these observables. Many other aspects are thinkable and interesting 
but have to be postponed to future studies. 



A. Hydrodynamical expansion and hadronization 



We first briefiy discuss the evolution and hadronization of the QGP-cylinder present on 
the T = Ti hypersurface as obtained from the hydrodynamical solution. Similar arguments 
and results can be found in a variety of papers, see for example |T8| , p!9|j26| , |i3|j46| , |i9|j50| , |55[] . 

In particular, ref. employed the very same model as here (i.e. longitudinal scaling 
fiow with cylindrically symmetric transverse expansion, the initial conditions and the EoS). 
However, the evolution at CERN-LHC energy had not been covered, and the hadronization 
hypersurface was only shown for a step-function like initial transverse energy density dis- 
tribution, but not for the wounded-nucleon distribution employed here. Therefore, a short 
discussion of the prehadronic stage may be in order here. 

Fig. D summarizes the space-time picture in the plane rj = (j) = 0. We show projections 
of various hypersurfaces on the (r, r7-)-plane because their shape in the (p- and //-directions 
is trivial: they are simply horizontal lines in the (r, (p)- (r, ■?7)-planes, extending from — vr to 
vr and —00 to 00, respectively. Thus, no derivatives like dr/dcp etc. appear in da^, eq. (pl|). 

Basically, we start at r = with a pure QGP extending from = up to tt = Rt = 
6 fm in the transverse direction. The thickness of the non-QGP region at the surface is very 
small for the wounded nucleon distribution. Initially, the hot quark-gluon fiuid is cooled 
mainly due to the longitudinal expansion, except close to the surface, where transverse 
pressure gradients are also large and lead to expansion rates several times larger than the 
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simple 1/r law [^,^. The fluid eventually reaches the boundary to the mixed phase, 
denoted by A = 1. (A is the local fraction of quarks and gluons within the mixed phase.) 
Clearly, the space-time volume of pure QGP increases substantially from SPS to RHIC and 
then again towards LHC This leads to stronger transverse flow of matter entering the mixed 
phase at RHIC and LHC than at SPS. Due to this effect the hadronization hypersurface 
(A = 0) extends to larger r^. 

At SPS, the hadronization hypersurface A = 0, where the switch to the microscopic 
model is performed, is almost stationary for some time r ^ Rt, after which the entire 
fluid hadronizes rapidly. UrQMD is being fed with hadrons from the stationary surface of a 
"burning log" of mixed phase matter ||50|. This is not to be misunderstood as an evaporation 
process, though. The fluid is moving with substantial velocity through the hadronization 
hypersurface, in particular near the point where the A = and the r = hypersurfaces 
meet (there, the very dilute fluid comes close to the light-cone). Thus, the momenta of the 
emitted hadrons, which are purely thermal in the local rest frame, are boosted in transverse 
direction. At RHIC, and of course even more so at LHC, that pre-acceleration by the 
QGP "explosion" is so strong that the hadronization hypersurface is initially even driven 
outwards, before the mixed-phase cylinder flnally collapses (when it can not balance the 
vacuum pressure B any more) and emits hadrons from all over the transverse area. 

Thus, it is clear from Fig. ^ that the dynamics at SPS is characterized by the large 
space-time volume occupied by the mixed phase, while the stiffer and more "explosive" 
|50,|73| QGP gains importance at higher energies. 



B. Post-hadronization kinetics: evolution of (px) 



The choice of the hypersurface at which to perform the transition from the macroscopic 
hydrodynamical calculation to the microscopic transport model may affect the reaction 
dynamics and the results of our calculation. However, concerning the variation of the hy- 
persurface for that transition one has to note that the hadronic part of the EoS used in 
the hydrodynamic solution contains the same states as UrQMD, and the energy-momentum 
tensors on both sides of the hadronization hypersurface match. // the assumption of local 
equilibrium is indeed fulfllled, UrQMD will simply continue the hydrodynamic flow since it 
reduces to hydrodynamics in the equilibrium- limit. However, as we shall see later, for some 
hadron species with small interaction cross sections deviations from ideal hydrodynamic flow 
can be observed immediately after complete hadronization (see also refs. [^,0]). It is found 
that the expansion of the hadronic fluid is dissipative rather than ideal; due to the fast local 
expansion generated by the QGP before hadronization the ideal flow is disturbed. Therefore 
it does not make much sense to choose a later hypersurface for the matching because one 
would precisely assume that ideal flow persists even after hadronization. 
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Nevertheless, it is interesting to study how the choice of a later hypersurface for the 
transition from the macroscopical to the microscopical part of the calculation affects the 
results. This reveals "how wrong" the assumption of an ideal evolution of the state at 
hadronization is. Figure ^ shows the final mean transverse momentum (pt) for various 
hadron species as a function of the temperature on the hydro— s-micro transition isotherm, 
Tsw The grey lines in the upper frame denote the {pr) of the hadrons at hadronization, i.e. 
at Tc = 160 MeV. As shall be discussed in greater detail in section |1V li| , the change in {pt) 
in the hadronic phase (for our "default" choice Tgy^ = Tc = 160 MeV) depends strongly on 
the individual hadron species. Protons and hyperons gain most, the Q~ does not acquire 
any additional (pt) at all, and pions even loose some (px) due to rescattering and additional 
soft pion production. 

The results change only marginally when decreasing T^^ to 150 MeV. Simply speaking, 
UrQMD reproduces the fluid-dynamical solution down to about T ^ 150 MeV, for central 
Au+Au at RHIC energy. At this stage, fluid-dynamics predicts that the transverse rarefac- 
tion in the hadron fluid reaches the center. Consequently, the expansion becomes rather 
spherical and transverse flow increases strongly in this "hadronic explosion". The lower 
frame in Fig. § shows the kink in {pt){Tsw) of heavy hadrons at Tgw ~ 150 MeV predicted 
by ideal hydrodynamics. 

Remarkably, however, the system apparently is already in a state of too rapid expansion 
for this "hadronic explosion" to happen. Given the state at hadronization, UrQMD (applying 
realistic cross-sections) predicts that the hadronic fluid basically freezes out right at the point 
where the hadronic rarefaction is about to make the expansion more spherical and to increase 
the expansion rate, see e.g. Fig. 2 in |^. Any later transition from hydrodynamics to the 
microscopic transport model leads to a strong increase of (pt) at freeze-out, which depends 
only on the mass of the hadron, but not on its flavor (resp. its quark content). 

The lines in the lower frame of figure § show the (px) of the respective hadron species at 
the transition hypersurface (i.e. at Tsw)- By comparing the (pt) value indicated by the line 
to that given by the plot symbol for each Tsw one can determine the amount of (px) gained 
or lost during the microscopic evolution of the reaction. Again, protons acquire the most 
{Pt) during the microscopic evolution (even though the amount of {pr) gained decreases 
the lower Tsw is and the closer the system comes to freeze-out), whereas S's and f2's do not 
experience any (px) increase at all. 

It is obvious from this analysis that the conditions of applicability for hydrodynamics 
in the hadronic phase deteriorate rapidly. A general freeze-out criterion can not be given 
since the freeze-out depends on the system size and the centrality, the energy etc. However, 
our transport calculation with realistic cross-sections in the hadron gas, starting in the 
wake of a hadronizing QGP, shows that the expansion is too rapid to allow cooling of the 
strong interactions much below Tc- In particular, adiabatic expansion breaks down once 
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the expansion of the hadron fluid effectively becomes (3+l)-dimensionaL 



C. Space-time distributions of hadronic freeze-out 

Let us now turn to the freeze-out "hypersurfaces" of pions and nucleons in central (impact 
parameter 6 = fm) collisions of gold or lead nuclei at SPS {^/s = 17 GeV per incident 
colliding nucleon-pair) , RHIC {^/s = 200 GeV per incident colliding nucleon-pair) and LHC 
{^/s = 5500 GeV per incident colliding nucleon-pair). We start with the nucleons, the most 
abundant baryon species in the system, restricting ourselves to the central rapidity region. 
Figure]^ shows the freeze-out]^ time and transverse radius distributions I/tt d^A^/drydrj^dy 
for LHC (top), RHIC (middle) and SPS (bottom). The right column shows the result of the 
pure hydrodynamical calculation up to complete hadronization, with subsequent hadronic 
resonance decays, but without hadronic reinteraction. The left column shows the same 
calculation including full microscopic hadronic dynamics. 

The freeze-out characteristics of the nucleons are signiflcantly modifled due to the 
hadronic interaction phase. The average transverse freeze-out radius doubles at SPS and 
RHIC and increases by a factor of 2.5 at LHC (see also table |T|). The respective aver- 
age freeze-out times increase by similar factors (see table [TV] ). E.g., at RHIC the average 
freeze-out time for protons changes from 11.3 to 25.8 fm/c due to hadronic rescattering. 

As the meson multiplicity in the system at RHIC is flfty times larger than the baryon 
multiplicity, baryons propagate through a relativistic meson gas, acting as probes of this 
highly excited meson medium. Thus, we use the proton and hyperon freeze-out values listed 



in table |^ for a flrst rough estimate of the duration of the hadronic phase via AT^ad = 

^^Hydro+UrQMD^ _ ^^Hydro+had.decays^ _ gpg ^^^^^ ^^^^^ ^ 5 ^^^^^ ^^^^ 

similar to the value at RHIC (~ 15 fm/c) and at the LHC we obtain Arhad ~ 23 fm/c. 
The transverse spatial extent of the hadronic phase can be estimated in a similar way, using 

Hydi 



table ^ and deflning the thickness Avhad of the hadronic phase as: Avhad = (^r^J^''°+U''Q^°^ _ 



('"f/r™"'^'^^'^ '^'^'^'^'^'^)- Here we flnd values of ~ 4.4 fm at the SPS, ~ 5.8 fm at RHIC and 
13.3 fm at the LHC. 

The Hydro-l-UrQMD model predicts a space-time freeze-out picture which is very differ- 



ent from that usually employed in the hydrodynamical model, e.g. in refs. ||2T| , |25| , P6| , |46| , |74 



Here ||3^, freeze-out is found to occur in a four- dimensional region within the forward light- 



cone 1 31 1 rather than on a three-dimensional "hypersurface" [38|. Similar results have also 



^Preeze-out meaning the space-time point of last interaction, irrespective of how "soft" that last 
interaction might be. We remind you also that mean-fields are not taken into account. They could 
even prolong the freeze-out due to very soft interactions of the hadrons with the mean field. 



23 



been obtained within other microscopic transport models when the initial state was 
not a quark-gluon plasma. This finding seems to be a generic feature of such models: the 
elementary binary hadron-hadron interactions smear out the sharp signals to be expected 
from simple hydro. This predicted additional fourth dimension of the freeze-out domain 
could affect the HBT parameters considerably. 

This does not mean that the momentum- distributions alone can not be calculated assum- 
ing freeze-out on some effective three-dimensional hypersurface. For example, if interactions 
on the outer side of that hypersurface are very "soft" , the single-particle momentum distribu- 
tions at not too small pt will not change anymore. The two-particle correlator does change, 
however, since it probes rather small relative momenta. Thus, the freeze-out condition, e.g. 



the temperature, as measured by single-particle spectra and two-particle correlations ||75 
needs not be the same. 

The shapes of the freeze-out hypersurfaces (FOHS) show broad radial maxima for inter- 
mediate freeze-out times. Thus, transverse expansion has not developed scaling-flow (in that 
case the FOHS would be hyperbolas in the t — tt plane). This agrees with the discussion 
of the evolution of the {px) after hadronization in section |1VB , which already indicated the 
transition to free-streaming once the transverse expansion rate becomes comparable to the 
longitudinal expansion rate. 

Furthermore, the hypersurfaces of pions and nucleons, and their shapes, are distinct from 



each other (as also found in |2|j3|j3|j4g] at the lower BNL-AGS and CERN-SPS energies). 
Thus, the ansatz of a unique freeze-out hypersurface for all hadrons appears to be a very 



rough approximation, cf. also refs. ||32| , |33| ,|3< 

Figure |^ shows the transverse freeze-out radius distributions for tt, ii', A + S and 
Vt^ at LHC (top), RHIC (middle) and SPS (bottom). They are rather broad and similar to 
each other, though the VL~ shows a somewhat narrower freeze-out distribution. The average 
transverse freeze-out radii are listed in table |T|; e.g. at RHIC we find 9.5 fm for pions, 
10.2 fm for kaons, 11.3 fm for protons, 11.6 fm for Lambda- and Sigma- Hyperons, 14.2 fm 
for Cascades, but only 7.3 fm for the Vt^ . The freeze-out of the VL~ occurs rather close 
to the phase-boundary due to its very small hadronic interaction cross section. This 
observation holds true for all three studied beam energies. The respective thickness Arhad 
of the hadronic phase is reduced by a factor of 2 for the Q~ , compared to that of the other 
baryon species. This behavior could be responsible for the experimentally observed hadron- 
mass dependence of the inverse slopes of the m-r-spectra at SPS energies p6|. For the fl~, 
the inverse slope remains practically unaffected by the purely hadronic stage of the reaction, 
due to its small interaction cross section, while the flow of p's and A's increases (see 
also section |LV hp . 

Figure |^ shows the freeze-out time distributions d'^N/drfrdy for vr, p and at LHC 
(top), RHIC (middle) and SPS (bottom). Open symbols denote the distributions for a pure 
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hydrodynamical calculation up to hadronization with subsequent hadron resonance decays 
(but without hadronic reinteraction) , whereas the full symbols show the full calculation 
with hadronic rescattering. As we have already seen previously in the transverse freeze-out 
radii, hadronic rescattering strongly modifies the shape of the distributions and significantly 
increases the lifetime of the system. Table |V| lists the average freeze-out times for tt, p, 
Y{= A + S°), S and Q~ with and without hadronic rescattering. 

One issue of great interest is the predicted significant increase of the lifetime of the 
system from SPS to RHIC energies |Q, being due to the time-delay caused by a first-order 
phase transition [^. However, our model calculation (which does exhibit a first order phase 
transition) shows no huge difference in the freeze-out time distributions of vr, p, and Q~ from 
SPS to RHIC energies (note, however, the logarithmic scale). Origin of this prediction is 
that we include many more states in the hadronic EoS, which speeds up hadronization 
considerably Furthermore, decays of resonances partly hide the remaining small 

increase of the hadronization time. Thus, the "time-delay signal" can not be expected to be 
well above ~ 20 — 30%, and must be approached by a detailed excitation function. 

Note that the multi-strange Q~ baryons freeze out far earlier than all other baryons, 
as discussed already previously in the context of figure p. The duration of the hadronic 

. 1 !• 1 A / Hydro+UrQMD\ / Hydro+had. decays \ • i i j 

remteraction phase, /AThad = [Tfr } — {Tj^ ) remams nearly unchanged, 

e.g. at 5.9 fm/c for pions, 8.0 fm/c for kaons, 14.5 fm/c for protons, 15.4 fm/c for hyperons 
and 8.0 fm/c for the VL~ between RHIC and SPS. 

Note that the lifetime of the pre-hadronic stage in this approach is a factor of 2 — 3 
longer than when employing the parton cascade model (PCM) ||77| , |78| for the initial reaction 
stage. It will be interesting to check whether this is related to the first-order phase transition 
built into the EoS which is used here. The final transverse freeze-out radii and times (after 



hadronic rescattering), however, are very similar in both approaches |[78 



Figure || shows the estimated freeze-out volume V* = 'n'{rTjr)^'Tfr as a function of the 
pion rapidity density dA^/dy for four different bins in transverse momentum. For all pt- 
bins V* exhibits a nearly linear increase with dN-^/dy. Thus, the freeze-out density of the 
pions remains virtually constant over a large range of multiplicities (or energies). We will 
see in the next section that this is due to the fact that the chemical freeze-out of pions 
occurs rather shortly after hadronization of the QGP, at all energies studied here. Since the 
local density of pions on the hadronization hypersurface is similar in all cases (because the 
temperature is almost the same), the density at chemical freeze-out is, too. 

We also observe that low-pr pions are basically emitted from the entire volume, while at 
higher p^ the pions seem only to be emitted from an outer shell, the radius of the hollow core 
increasing with pt- The inset of figure ^ shows the dependence of the non-emitting core- 
volume Vo on the transverse momentum of the pions. Vq has been calculated by a linear fit 
of V* to dA'Tr/d?/: V^* = Vq + c{dNT,/dy). The increase of Vq with px is a manifestation of 
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the collective flow effect; high-pr pions can not be emitted from the center, ~ 0, since 
the collective velocity fleld vanishes there. 



D. Chemical freeze-out 

So far, we have only discussed the kinetic freeze-out of individual hadron species. How- 
ever, apart from the kinetic freeze-out, the chemical freeze-out of the system, which fixes 
the chemical composition, is of interest, too. 

The chemical freeze-out hypersurface of hadron species i is in principle defined as the 
surface o"^/jem separating the space-time region where d-Ni = from that where the number- 
current Ni is not conserved. Usually, the chemical freeze-out is defined modulo hadronic 
resonance decays which are performed on c^/jg^, even for short-lived resonances like the p- 
meson or A-baryon. However, that definition is not very useful in the present case, since 
most inelastic processes are actually modeled via resonance excitation and subsequent decay, 
cf. section [HI D| . Furthermore, as in the case of kinetic freeze-out studied above, the micro- 



scopic transport model does not yield sharp hypersurfaces (three-dimensional volumes) but 
rather freeze-out domains (four-dimensional volumes). We shall therefore mainly discuss the 
evolution of hadron multiplicities after hadronization, and their time-dependence. 

Figure ^ shows the time evolution of on-shell hadron multiplicities for LHC (top), RHIC 
(middle) and SPS (bottom). The dark grey shaded area indicates the duration of the QGP 
phase, whereas the light grey shaded area depicts the mixed phase (both averaged over r^; 
only hadrons that have already "escaped" from the mixed phase into the purely hadronic 
phase are shown). Hadronic resonances are formed and are populated for a long time. 
One can rather nicely observe the stronger transverse expansion as beam energy increases: 
on r = const, hypersurfaces the resonance- decay "tails" get boosted to larger r. Due to 
those transversely boosted resonances the hadron yields saturate only at rather large r, 
approximately 25 fm/c at SPS and RHIC and about 40 fm/c at LHC. 

By comparing the final hadron yields resulting from the hydrodynamical calculation (up 
to hadronization, including subsequent hadronic decays, but no hadronic reinteractions) to 
that of the full calculation, which includes microscopic hadronic dynamics, we can quantify 
the changes of the hadrochemical content due to hadronic rescattering. 

Figure |10| shows the relative change (in %) of the multiplicity for various hadron species 
for SPS (bottom), RHIC (middle) and LHC (top). As to be expected, the state of rapid 
expansion prevailing at hadronization does not allow chemical equilibrium to hold down 
to much lower temperatures. The hadronic rescattering changes the multiplicities by less 
than a factor of two, cf. also [|1^]. Thus, we have first evidence that a QGP expanding and 
hadronizing as an ideal fiuid produces a too rapidly expanding background for a hadron-fiuid 
with known elementary cross-sections to maintain chemical equilibrium down to much lower 
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temperatures than Tc- 

However, a closer look provides more insight into the chemical composition. The changes 
are most pronounced at the SPS, were the baryon-antibaryon asymmetry is highest (since 
the net-baryon density at mid-rapidity is highest). This manifests e.g. in a reduction of 
the antiproton multiplicity by 40-50% due to baryon-antibaryon annihilation. A and S are 
affected in similar fashion. 

The baryon-antibaryon asymmetry decreases at higher beam energy, and at LHC particle- 
antiparticle symmetry is restored for our initial conditions. The remaining small asymmetries 



(compare e.g. the p-p, K-K, and Y-Y evolutions in Fig. [1^) are due to fluctuations triggered 
by the finite number of particles, which distort the ideal longitudinal boost-invariance present 
(by construction) at hadronization. 

Interestingly, the VL~ multiplicity decreases stronger towards higher beam energy. This 
is due to the higher antibaryon density in the system, leading to more VL~ annihilations on 
antibaryons with subsequent redistribution of the three strange quarks. (This process is 
modeled in UrQMD as string excitation and subsequent fragmentation, cf. |^.) Thus the 
hadronic phase becomes slightly more opaque for the Vt~ with increasing beam-energy. 

Collision rates offer another approach to determine the duration of the hadronic phase, 
in particular B-B collisions which almost always lead to annihilation. Fig. ^ shows the 
time-evolution of the rates for hadron-hadron collisions at LHC (top), RHIC (middle) and 
SPS (bottom). Meson- meson (MM) and - to a lesser extent - meson-baryon (MB) in- 
teractions dominate the dynamics in the hadronic phase at RHIC and LHC. At the SPS 
meson-baryon and meson- meson interactions are equally frequent. Note that while at the 
SPS baryon-baryon (BB) collisions significantly outnumber baryon-antibaryon annihilations, 
the situation at RHIC and LHC is reversed, where B-B annihilation is far more frequent 
than BB collisions. This is a consequence of the fact that the B-B annihilation cross sections 
at small relative momenta increase faster then the total B-B cross sections [^. In the case 



of (approximate) baryon-antibaryon symmetry, one therefore expects more B-B than B-B 
interactions, as seen for RHIC and LHC energies. 

Of course, all collision rates reach their maxima at the end of the mixed phase, then 
decreasing roughly according to a power-law. After 35 fm/c, less than one hadron-hadron 
collision occurs per unit of time and rapidity at SPS and RHIC energies; due to the higher 
transverse 7-factor the time is ~ 60 fm/c at the LHC. At this stage the system is certainly 
kinetically and chemically frozen-out. 

E. Transverse flow: Emission of multi-strange baryons from the phase-boundary 

In this section we analyze the transverse mass spectra at freeze-out, and discuss their 
evolution from the hadronization hypersurface. The results obtained for Pb+Pb collisions 
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at CERN-SPS energy are in reasonable agreement with the data obtained by the NA49- 



collaboration and by the WA97-collaboration [^]. For a comparison to that data we 
refer to here, we focus on the model-results. 

Fig. |12| compares the m^^-spectra on the hadronization hypersurface (open symbols), ob- 
tained from Eq. ( ^21) (plus strong resonance decays), with those at freeze-out (full symbols). 
One observes that the transverse flow of p's and A's increases during the hadronic stage, since 
those spectra flatten. On the other hand, the spectra of f2's and of S's with mrp > 1.6 GeV 
are practically unaffected by the hadronic stage and closely resemble those on the phase 
boundary. This is due to the fact that the scattering rates of S and Q in a. pion-rich hadron 
gas are significantly smaller than those of A^'s and A's ||36|j37|j80|| . As shown in Fig. on 
average the baryons which finally emerge as S's and il's suffer far less interactions than the 
final-state p's and A's. Thus, within the model presented here, these particles are basically 
emitted directly from the phase boundary with very little further rescattering in the hadronic 
stage. The hadron gas emerging from the hadronization of the QGP (in these high-energy 
reactions) is almost "transparent" for the multiple strange baryons. On the other hand, p's 
and A's on average suffer several collisions with other hadrons before they freeze-out. This 
behavior holds generally true for all three studied energy domains, at the SPS, RHIC and 
LHC. 

These findings manifest themselves most strikingly in the mass-dependence of the inverse 
slopes of the rriT spectra. A simple isentropic hydrodynamical expansion leads to broader 
rriT spectra of heavier states, i.e. (pt) or the inverse slope T* increase with mass |PT|. This 
observation agrees with the inverse slopes of vr, K, and p measured for central collisions of 
Pb nuclei at a CM-energy of 17 A GeV |^. However, it has also been found that the S and 
Q baryons do not follow this general trend [0,^ . 

Fig. depicts the inverse slopes T* obtained from our model by a fit of d'^Ni/d'^rriTdy 
to exp(— m^/T*) in the range — rrii < 1 GeV. The statistical error of this fit is ~ 10%. 
Open symbols denote the SPS calculation and data, whereas full symbols show the RHIC 
prediction. The lines show a purely hydrodynamical calculation pT| , |37[| with a freeze-out 
temperature of Tjo = 130 MeV for SPS (dotted line) and RHIC (full line), respectively. The 
trend of the SPS data (open circles), namely the "softer" spectra of S's and fi's as compared 
to a linear T*(m) relation, is reproduced reasonably well. As already mentioned, this is not 
the case for "pure" hydrodynamics with kinetic freeze-out on a common hypersurface (e.g. 
the T = 130 MeV isotherm), where the stiffness of the spectra increases monotonically 



with mass, cf. Fig. n4 and also refs. [26,54|. Resonance decays are not included in the 



hydrodynamic spectra on the T = 130 MeV isotherm. 

When going from SPS to RHIC energy, the model discussed here generally yields only a 
slight increase of the inverse slopes, although the specific entropy is larger by a factor of 4-5 ! 
The reason for this behavior is the first-order phase transition that softens the transverse 
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expansion considerably [0]. For our set of initial conditions, the average collective transverse 
flow velocity (at mid-rapidity) on the hadronization hypersurface increases only from ^ 0.3 
(for Pb+Pb at SPS) to ^ 0.35 (for Au+Au at RHIC) 0. (However, there are high- 
vt tails on the hadronization hypersurface which get more pronounced at RHIC than at 
SPS.) As can be seen from the present calculation, this is not counterbalanced by increased 
rescattering in the purely hadronic stage - compare to the inverse slopes obtained from 
"pure" hydrodynamics with freeze-out on the T = 130 MeV isotherm! 

The transverse flow at LHC beam energy is so strong that the rriT-spectra can not be 
fitted any more by an exponential distribution. We have therefore refrained to extract the 
slopes for the LHC calculation. Instead, in figure |1^ we show the mean transverse momenta 
of the different hadron species as a function of their mass. As in figure we compare the 
(Pt) on the hadronization hypersurface (open symbols), obtained from Eq. ( P^D (plus strong 
resonance decays), with that at freeze-out (full symbols). Hadronic rescattering leads to a 
transfer of transverse energy/momentum from pions to heavier hadrons (the pions actually 
suffer a reduction of (pt) in the hadronic phase) [^. This phenomenon has also been 
termed the pion-wind p6| , p3| , pushing heavier hadrons to higher pt- Nucleons gain most 
transverse momentum, while the Q~ remains nearly unchanged due to its small interaction 
cross section in the meson dominated hadronic medium, as discussed earlier in this section. 
Those hadrons are the best "messengers" of the early pre-hadronization evolution. 

Furthermore, one clearly observes the rather moderate increase of (px) from SPS to RHIC 
energy, as discussed already in Fig. II. In contrast, in our model the collective dynamics at 
the much higher CERN-LHC energy is dominated by the stiff QGP, cf. also Fig. |^, and the 
average transverse momenta increase appreciably. 



V. SUMMARY AND OUTLOOK 

In summary, we have introduced a combined macroscopic/microscopic transport ap- 
proach, combining relativistic hydrodynamics for the early deconfined stage of the reac- 
tion and the hadronization process with a microscopic non-equilibrium model for the later 
hadronic stage at which the hydrodynamic equilibrium assumptions are not valid anymore. 
Within this approach we have self-consistently calculated the freeze-out of the hadronic sys- 
tem, accounting for the collective flow on the hadronization hypersurface generated by the 
QGP expansion. 

The reaction dynamics, hadronic freeze-out and transverse flow in ultra-relativistic heavy 
ion collisions at SPS, RHIC and LHC have been discussed in detail. We find that the 
space-time domains of the freeze-out for the investigated hadron species are actually four- 
dimensional, and differ drastically between the individual hadrons species. 

The thickness of the hadronic phase is found to be between 2 fm and 6 fm (at RHIC), 
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depending on the respective hadron species. Its lifetime is between 5 fm/c and 13 fm/c, 
respectively. Freeze-out radii distributions have similar widths for most hadron species, 
though the Q~ is found to be emitted rather close to the phase boundary and shows the 
smallest freeze-out radii and times among all baryon species. The total lifetime of the system 
does not increase drastically when going from SPS to RHIC energies. 

Our model-calculation shows that in high-energy nuclear collisions the hadron multiplic- 
ities at midrapidity change by less than 40% after hadronization, unlike e.g. in the early 
universe. However, a closer look is warranted and reveals interesting information. For ex- 
ample, more strange baryons (A, S, S, Q) are annihilated as the energy increases because 
the anti-baryon density at hadronization increases. 

Interactions within the hadron gas increase the collective flow beyond that present at 
hadronization, and reduce the temperature below the QCD phase transition temperature (we 
assume Tc = 160 MeV). As an exception, we find that multiple strange baryons practically 
do not rescatter within the hadron gas. Their mr-spectra are therefore determined by the 
conditions on the hadronization hypersurface, i.e. Tc and the collective flow created by the 
expansion preceding hadronization. Their spectra therefore are less sensitive to the confined 
phase, T < Tc, but are closely related to the EoS of the QGP and the phase transition 
temperature Tc- 

Average transverse momenta and inverse slopes are predicted to increase only moderately 
from SPS to RHIC, despite the significant increase of the entropy to net baryon ratio. In this 
sense, the collective evolution at RHIC energy is strongly characterized by the presence of a 
well-mixed coexistence phase with small isentropic speed of sound. It will be very interesting 
to see if this picture of hadronization of bulk QCD matter, which is based on similar models 
for the QCD phase transition in the much slower expanding early universe, agrees with the 
data to be taken by the various experiments at BNL-RHIC. 

Towards the much higher CERN-LHC energy, the evolution changes appreciably. The 
pure QGP occupies a larger space-time volume than the mixed phase. If the QGP-EoS 
at high energy density is anywhere close to an ultrarelativistic ideal gas with p ~ e/3, 
transverse expansion should be much stronger than at RHIC and SPS. Consequently, the 
average transverse momenta of the heavier hadrons increase by 60 — 80% as compared to 
RHIC energy. 

We believe that the model presented here and in refs. represents a step forward 

towards the understanding and the description of the evolution of a quark-gluon plasma, 
its hadronization, and the subsequent freeze-out of the strong interactions. Nevertheless, 
it is clear that many improvements are thinkable and necessary before a really detailed 
comparison to experimental data can be attempted. 

For example, corrections to ideal fluid dynamics (before hadronization) should be studied, 
at least within the Navier-Stokes approximation. In the present approach dissipative effects 
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are only taken into account after hadronization, where we expect them to be most significant, 
particularly as freeze-out is approached. 

The widely used Bag-model EoS can certainly be improved as well. It is well known that 
it yields a substantially higher latent heat than extracted from present lattice-QCD results. 
Thus, it may over-pronounce the effects of a first-order QCD phase transition. One may 
even try a cross-over transition to see whether that is ruled out by experimental data or 
not. Also, we have already commented on the fact that due to the huge expansion rate in 
high-energy collisions (which is not much smaller than strong interaction rates) more radical 
scenarios like spinodal decomposition rather than an adiabatic phase transition should be 
examined as well. 

To simplify the switch from hydrodynamics to the microscopic transport model we as- 
sumed longitudinal boost-invariance and azimuthal symmetry. The latter approximation, 
in particular, disables us to study many up-to-date topics that will be addressed by the 
experimental data, as e.g. anisotropies in the hadron spectra (in non head-on collisions) 
that may be sensitive to the QGP-EoS. A fully 3+1 dimensional solution without symmetry 
assumptions is thus highly desirable. 

One application not covered at all in the present studies is correlated particle emission. 
Two-particle correlations may for example allow to extract the volume of the emission region 
in space-time. The two-particle correlator probes rather soft interactions and therefore 
contains information about the freeze-out process (e.g. the thickness of the emission region, 
background mean-fields etc.) ||75| , [75| . 

Fluctuations in the rapidity and transverse momentum spectra induced (or suppressed) 
by the QCD phase transition are another highly interesting topic. So far, we have studied 
only the evolution on average, and have allowed only for fiuctuations to develop in the post- 
hadronization stage (which occur naturally in the microscopic transport model). However, 
the hadronization process could in principle induce larger or other types of fiuctuations, e.g. 
due to droplet formation [Q, spinodal decomposition with possible DCC formation [^, or 
due to the change of the order of the phase transition in the vicinity of a second order critical 
point [^. It would obviously be highly desirable to know if they can survive the hadronic 
interaction stage. Other fiuctuations may be there from the very beginning, e.g. those 
arising in the production process of secondary hadrons, and one could study the evolution of 



strangeness-rich rapidity bins |£3[, or of rapidity bins with negative baryon- number, through 
the hadronization phase transition until freeze-out. 

There are many other questions that can not be listed here but can be addressed within 
this model. Work along those lines is in progress and will be reported in forthcoming 
publications. 
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TABLE I. Baryons and baryon-resonances treated in the model. The corresponding antibaryon 
states are included as well. 



0- 1- 0+ 1+ 2+ (1-)* 

TT p ao ai 02 /9(1450) 

K K* Kl p(1700) 

rj u h h /2 a;(1420) 

ri' <t> /o* /( /2 ^(1600) 



TABLE IL Mesons and meson-resonances, sorted with respect to spin and parity, treated in 
the model. 
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TABLE III. Mean transverse freeze-out radii (r^) for different hadron species at SPS, RHIC 
and LHC. d.p.b denotes the transverse freeze-out radii of the hydrodynaniical calculation up to 
hadronization, including subsequent hadronic decays, but no hadronic reinteractions, h.r. denotes 
the full Hydro+UrQMD calculation, including hadronic rescattering and h.t. stands for an estimate 



of the thickness Arhad of the hadronic phase: Arhad = ^7-^^'^''°+^'"'^'^^^ _ ^j^Hydro+had.decays^^ 
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TABLE IV. Mean freeze-out times for different hadron species at SPS, RHIC and LHC. d.p.b 
denotes the freeze-out times of the hydrodynamical calculation up to hadronization, including 
subsequent hadronic decays, but no hadronic reinteractions, h.r. denotes the full Hydro-l-UrQMD 
calculation, including hadronic rescattering and h. d. stands for an estimate of the duration of the 

hadronic reinteraction phase, Ar^ad = (Hydro+UrQMD^ _ ^ Hydro+had.decays^ 
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FIGURES 



time 




space 

FIG. 1. Schematic overview of the space-time evolution of a high-energy heavy-ion colUsion as 
assumed in the model presented here. 
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FIG. 3. Hypersurfaces corresponding to A = 1 (boundary between pure QGP and mixed 
phase) and A = (boundary between mixed phase and pure hadron phase) for LHC (top), RHIC 
(middle) and SPS (bottom). 
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FIG. 4. Mean transverse momentum {pt) of various hadron species at freeze-out (symbols) 

vs. the macro to micro transition temperature Tgw The horizontal lines in the upper frame show 

the respective (pt) values right after hadronization (A = hypersurface). The lines in the lower 

frame show the (px) emerging from ideal flow down to T = Tgw This figure is for central Au + Au 

collisions at BNL-RHIC energy. 
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FIG. 5. Freeze-out time and transverse radius distribution d^Ar/(rTdrrdr/rdy) for nucleons 
at LHC (top), RHIC (middle) and SPS (bottom). The left column shows the result for the pure 
hydrodynamical calculation up to hadronization with subsequent hadron resonance decays (but 
without hadronic reinteraction). The right column shows the analogous calculation, but with full 
microscopic hadronic collision dynamics after the hadronization. The contour lines have identical 
binning for all rows and columns. 
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at LHC (top), RHIC (middle) and SPS (bottom). 



46 



10* 



10" 



>%lt"' 



•a 



Iff' 
lo' 

10° 



>,iff' 

if 

5iff^ 



Iff' 



Iff* 



Iff' 

10 



10" 



Iff 



e 








01 



[5 [] 



fiill symbols: Hydro + UrQMD 

open symbok: Hydro -|- had. decay afrer hadrnnlzation 



RHIC 



• p 

■ V/5 

A S/IG 

♦ n/io 




H 1 1 1 1 1 1- 



5 



SPS 



• p 

■ Y/5 

A B/10 

♦ n/10 




** • 



10 20 30 to 50 60 70 
Tf, (fm/c) 



BO 



FIG. 7. Freeze-out time distributions d'^N/dTfj.dy of tt, p and il" for LHC (top), RHIC 
(middle) and SPS (bottom). Open symbols denote the distributions for a pure hydro dynamical 
calculation up to hadronization with subsequent hadron resonance decays (but without hadronic 
reinteraction) , whereas the full symbols show the full calculation with hadronic rescattering. 
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FIG. 8. Estimated freeze-out volume of pions as a function of the pion rapidity density for four 
different bins in transverse momentum. High p^-pions are only emitted from an outer shelly the 
radius of which increases with pr- The inset shows the volume of the "hollow core" as a function 
of PT- 
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FIG. 9. Time evolution of on-shell hadron multiplicities (integrated over rx) at LHC (top), 
RHIC (middle) and SPS (bottom). The dark grey shaded area shows the duration of the QGP 
phase whereas the light grey shaded area depicts the coexistence phase. 
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FIG. 10. Multiplicity change in % due to hadronic rescattering for various hadron species at 
SPS (bottom), RHIC (middle) and LHC (top). The error-bars give an estimate of the systematic 



error. 
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FIG. 11. Hadron-hadron collision rates at LHC (top), RHIC (middle) and SPS (bottom). The 
dark grey shaded area shows the duration of the QGP phase whereas the light grey shaded area 
depicts the coexistence phase. 
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FIG. 12. Transverse mass spectra of vr, K (left eolumn) and N, A+Tj^, H'^ + H", and 17^ (right 
column) at LHC (top), RHIC (middle) and SPS (bottom). The open symbols denote the spectra 
on the hadronization hypersurface whereas the full symbols show the calculation at freeze-out after 
hadronic rescattering. 
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FIG. 13. Distribution of the number of interactions that the final-state particles suffer after 
being hadronized; for LHC (top), RHIC (middle) and SPS (bottom). 
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FIG. 14. Inverse slopes of the mT-spectra of tt, K, p, k + TP, S'^ + and at ycm. = 0, 
rriT — rrii < 1 GeV. 
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FIG. 15. Average transverse momentum (px) vs. mass for tt, K, p, Y, S and $7. The open 
symbols denote the spectra on the hadronization hypersurface (including strong resonance decays), 
the full symbols show the value at freeze-out (after hadronic rescattering) . For clarity, the symbols 
have been shifted by ±10 MeV. 
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